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ABSTRACT 

A fully coherent acceleration search algorithm has been applied to the Parkes multi- 
beam pulsar survey of the Galactic plane to search for previously undiscovered rela- 
tivistic binary pulsars. The search has resulted in the discovery of 16 pulsars including 
a binary millisecond pulsar and an intermittent pulsar. Despite a number of promising 
candidates there have been no new discoveries of relativistic binary pulsars. Here we 
detail the acceleration search performed in our analysis and present coherent timing 
solutions for each of pulsars discovered. In light of the lack of discoveries of rela- 
tivistic binary pulsars, we also discuss the technique of acceleration searching and its 
effectiveness in finding these systems. 
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1 INTRODUCTION 

Radio pulsars in compact, highly-relativistic binary sys- 
tems, such as double neutron star systems (DNS) and 
pulsar white dwarf systems (NS-WD), provide the most 
precise tests of Einstein's theory of General Relativity 
(GR) and al ternative theories of gravi ty in the stron g-field 



regime (e.g. iTavlor fc Weisberd 1989 



ty m tlie strong-ncld 
Kramer et ai] 2006, 



iBreton et ail 2008, iFreire et all 2012). Future more strin- 
gent tests of gravity could be performed with more compact 
and/or higher-mas s systems such as pulsar- black hole binary 
systems (NS-BH) (|Kramer et all 2004), or even pulsars or- 
biting closely around the s uper massiv e blac k hole at th e cen- 
tre of the Galaxy, Sgr A* (|Pfahl et all 20Q5. lLiu et~ail 2012). 
The increased orbital velocities and deeper gravitational 
potentials present in such systems would result in larger 
post-Keplerian (PK) orb ital effects, which are meas urable 
through pulsar timing l|Damour fc Esposito-Faresel 1998, 
IWex fc Kopeikinl 1999). 

Searches for these highly-prized systems are notoriously 
difficult. After correction for the frequency dependent dis- 
persion in the arrival time of pulses, caused by the unknown 
free electron content along the line-of-sight ('dedispersion'), 
standard pulsar searches look for significant features in the 
fluctuatio n spectrum of the d edispersed time series (see, for 
example, iLorimer fc Kramerl 2005, for a detailed descrip- 
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tion of pulsar search methods). In relativistic binary sys- 
tems, Doppler smearing of the pulse frequency, caused by 
rapid orbital motion during the survey observation, renders 
standard Fourier-based pulsar search algorithms ineffective 
(see, e.g. Ijohnston fc Kulkarnil 1991). The binary search al- 
gorithms that compensate for these effects are computation- 
ally expensive in blind pulsar searches. Hence, the analysis of 
large-scale pulsar surveys, like the Parkes multi-beam pulsar 
survey (PMPS), have not yet be en performed with optimum 
binary searches (| Faulkner et all 2004) . In this work we have 
attempted to further combat binary selection effects in the 
Parkes multi-beam pulsar survey by re-processing the entire 
survey with an efficient 'coherent acceleration search'. 



The outline of the rest of this paper is as follows: In 
Section 2, we give a summary of the PMPS and describe 
previous binary searches of the survey data. Section 3 is a 
description of the binary search algorithm employed in this 
search. In Section 4, details of the re-processing data pipeline 
are given. Section 5 presents timing solutions of all 16 pulsars 
and further details on any of the notable discoveries. Section 
6 provides a discussion of possible reasons for the lack of 
discoveries of any new relativistic binaries in addition to 
further tests of the acceleration algorithm presented here. 
Section 7 is a summary and a discussion of future work on 
searches for relativistic binary pulsars. 



2 THE PARKES MULTI-BEAM PULSAR 
SURVEY 

The PMPS is the most successful pulsar survey ever com- 
pleted, havin g discovered 770 pulsar s up to the start 
of this work jManchester et all 2001 iMorris et ail 2002, 
Kramer et ail 2003, iHobbs et ail 2004, IFaulkner etall 2004, 



Lorimer et ail 2006, iKeith et ail 2 009) and over 30 p ulsars 



discovered b y thei r single puls es (|McLaughlin et ail 2006, 
iKeane et ail 2010, iKeane et ail 2011). The survey was car- 
ried out between 1997 and 2003 using a 13-beam receiver 
with a bandwidth of 288 MHz centred on 1374 MHz and 
deployed on the 64-m Parkes radio telescope in Australia. 
It covers a strip along the Galactic plane from I = 260° to 
I — 50° and |6| < 5° with time series of 35 minutes digitally 
sa mpled every 250 fi s. Full details of the survey can be found 
in lManchester et ail (2001). 



2.1 Previous binary searches 

Despite the early successes of the PMPS the number of com- 
pact binary pulsars (orbital period, Pb < 1 day) and mil- 
lisecond pulsars (spin period, P < 30 ms) discovere d was 
low; one and three respectively (Table 3.2. IFaulkner! 2004). 
The long integration time of 35 minutes gave the PMPS 
unprecedented sensitivity to long-period pulsars, although, 
in searches for relativistic binary pulsars, where the spin 
frequency can change over the duration of the observation, 
sensitivity was lost in standard Fourier periodicity searches 
due to smearing of the spectral features of the pulsar. 

After an improvement in computational resources at 
the Jodrell Bank Observatory, in the form of the CO- 
BRA 182-node Beowulf cluster, a full re-processing of the 
PMPS was performed using two search algorithms designed 
to address the previously-descr ibed binary selection effects 
(|Faulknerl 2004, IFaulkner et~ail 2004). The algorithms used 
in t his analysis w ere the 'stack-slide search' (see, for exam- 
ple^ Woo^^_ai| 1991) and the 'phase modulation search' 
l|Ransom et ail 2003). In the stack-slide search technique 
each dedispersed time series is split into a number of con- 
tiguous segments (16 in the case of the PMPS re-processing), 
which are Fourier transformed separately. The resulting fluc- 
tuation spectra are then summed together with different 
linear offsets to correct for linear drift in pulse frequency 
caused by any constant line-of-sight acceleration. So called 
'acceleration searches' have been shown to be valid for ob- 
servations whe re the integration time, T b s , is < ft,/ 10 
|Ransom et ail 2003). For the PMPS this corresponded to 
binary systems with Pb > 6 hr. 

The phase modulation search works by detecting "pe- 
riodic sidebands" in the fluctuation s pectrum; the co mmon 
marker of a pulsar in a binary system jRansom et ail 2003). 
Pulsars in binary systems where at least 1.5 orbits are com- 
pleted during the observation time can be detected with this 
technique. For the PMPS this corresponded to systems with 
P b < 0.4 hr. 

The IFaulkner et~ail (2004) reprocessing found 11 bi- 
nary pulsars including the double neutron star system, PSR 
J1756— 2251, that wo uld not have be en discovered without 
the stack-slide search (|Faulkner et ail 2005). Of the other 10 
binary systems, three were detected with higher significance 
because of the stack-slide search (on average by a factor of 



~ 1.9: Table 7.2, iFaulknerl 2004). Results from the phase 
modulation search have never been fully inspected due to 
difficulties in reconstructing the predicted orbit of pulsar 
candidates (Faulkner, private communication). 

Although the stack-slide search method improved sensi- 
tivity to pulsars in some fast binary systems, sensitivity was 
lost due to the incoherent summation of spectra. By split- 
ting the time series into a number of segments, each of which 
is Fourier transformed and then summed, the phase coher- 
ence of any real signal is not preserved throughout the entire 
length of the observation. In the presence of Gaussian noise, 
and assuming stacks of equal length, the signal-to-noise ra- 
tio (S/N) of a real signal increases oc \fn wh ere n is the 
number of spectra stacked (|Bradv fc Creightonl 2000). As al- 
ready noted, the stack-slide and phase modulation searches 
left a gap in sensitivity for pulsars in binary systems with Pb 
between ~ 0.4 and ~ 6 hr. For systems that lie in this region 
of parameter space, sensitivity might also have been lost be- 
cause of observed changes in acceleration ('jerk') that the 
linear offsets applied in the stack-slide search cannot com- 
pensate for. Both the orbital-period sensitivity gap and the 
reduction in sensitivity due to the incoherent nature of the 
stack search might have been critical in the detection of a 
relativistic binary pulsar that was near the detection thresh- 
old of the survey. In this work we describe further efforts to 
tackle binary selection effects in the PMPS. 



3 A COHERENT ACCELERATION SEARCH 
ALGORITHM, PMACCN 

It is possible to perform acceleration searches in the time 
domain that are fully coherent and without the loss of sen- 
sitivity encountered in stacking algorithms. One such tech- 
nique is 'time domain resampling'; a method that enabled 
the discovery of nine binary millisecond pulsars in observa- 
tions of 47 Tucanae at 20-cm (ICamilo et ail|2000l ). Following 
the description in Section 2 of lCamilo et ail (2000), the time 
series is transformed into a frame which is inertial with re- 
spect to the pulsar by application of the Doppler formula 
to relate a time interval in the pulsar frame r to the time 
interval in the observed frame t. 



'{t) = TO (1 + V(t)/ C ) , 



(1) 



where v(t) is the changing line-of-sight velocity of the pul- 
sar, c is the speed of light, and the constant to is cho- 
sen such that r is equivalent to the sampling interval, 
Tsamp at the midpoint of the observation. New samples 
are computed from a linear interpolation over the origi- 
nal time series (|Middleditch fc Kristianlll984l ). Following re- 
sampling, the time series can be searched with standard 
Fourier techniques, as if th e pulsar were isolated (see for 
example lLorimer fc Kramerl 2005). To fully remove the ef- 
fects of orbital motion, the exact form of v(t) must be 
known, although in a blind search, where the orbital ele- 
ments are a priori unknown, this would require a search 
over five Keplerian parameters; a task that would re- 
quire significant compu tational resources (see for example, 
iDhurandhar fc Vecchid 2001, and discussion in Appendix 
A). To reduce the amount of computation required, the as- 
sumption of a constant line-of-sight acceleration, a, can once 
again be applied, so that v(t) — at. Now, searches in just 



one parameter, constant line-of-sight acceleration, a, can be 
made, but in a fully coherent manner. 

A step size between each acceleration trial, 5a was deter- 
mined by fixing the amount of signal advance or retardation 
with respect to either end of the integratiorQ to a value that 
was deemed acceptable. This value, which we shall now term 
'pulse broadening', r acC n, was chosen to be four times the 
sampling interval, r samp , i.e. if the pulsar signal fell between 
two acceleration trials this would be the maximum value of 
pulse broadening at the end points of the data. Because of 
the quadratic nature of the pulse broadening with time, due 
to a constant accelera tion, a (see for example Equation 4 in 
IJohnston fc Kulkarnil 1991), 

Taccn(t) = at* /2 C , (2) 

at least 50 per cent of pulses in the integration would be 
smeared in time by less than one time sample. Setting the 
pulse broadening time, given by Equation 2, equal to 8 r sam p 
allows a maximum of 4 r samp pulse smearing for signals that 
lie exactly in-between acceleration trials. Remembering that 
resampling is performed with respect to the midpoint of the 
data, t = T bs/2, giving the acceleration step size, 5a: 

5a = 64cT sa mp/T obs . (3) 

Our choice of acceleration step size is illustrated in Figure [T] 
Here we plot the pulse phase as a function of time for an 
idealized top hat pulse that has a width equal to a single 
sample, and that is accelerated to a value that is directly 
in between acceleration trials. The pulse phase broadening 
at the end points of the integration is therefore 4T ssmp ; the 
maximum value it can take. For a broadening of 4r samp at 
the end points it can be seen that 50 per cent of the pulses 
are smeared in phase by less than r aamp . Although the total 
broadening time is less than T aamp for half of the observation, 
the majority of power in a sample will be smeared into the 
next sample when, for example, r accn > 0.5r samp , which 
occurs after T_b_/\/32 using this step size. 

The acceleration step size given in Equation[3]scales lin- 
early with the sampling interval and with inverse square of 
the integration time. Longer integrations and smaller sam- 
ple intervals therefore require finer acceleration steps. Both 
of these parameters could be tuned within the limits of the 
PMPS data to find an optimum acceleration step size which 
covered new regions of parameter space without excessive 
levels of computation. The final processing parameters of 
our search algorithm, termed PMACCN, are given in Ta- 
ble [T] Data processing benchmar ks and algorithm perfor- 
mance tests described in lEatoughl (20090 favoured a search 
based on two independent half segments of the original 35 
minute integration, and a sampling interval increased by a 
factor of four from the original 250 /is to 1 ms. Using this 
choice of processing parameters a wide range of accelerations 
(±500 m s -2 ) could be searched in a manageable time on the 
computational hardware available (see Section 4). The ac- 
celeration range chosen is nearly twice the maximum accel- 
eration present in the most relativistic b inary pulsar kn own: 
the Double Pulsar system (~ 260 m s~ 2 : lBurgav et all 2003, 

1 The effect of the constant, to in Equation 1 is such that resam- 
pling is performed with respect to the midpoint of the observation. 

2 www.jb.man.ac.uk/~reatough/reatough_thesis.pdf 
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Figure 1. Time versus pulse phase for an idealized top hat pulse 
with a width equal to the sampling interval (dark grey blocks) 
and a true acceleration, indicated by the thick black line, that 
lies directly in between acceleration trials given by Equation [3] 
(thin black lines). Arbitrary time intervals are marked by the 
vertical dotted lines. The size of the acceleration step is indicated 
by the vertical arrow on the right. For signals that have a true 
acceleration that lies in between acceleration trials, at least 50 
per cent of the pulses are smeared by less than T samp . 



Table 1. PMACCN processing parameters. The step size in dis- 
persion measure ( DM) was determined following standard pro- 
cedures outlined in lLorimer fc Kramer] (2005). The range in DM 
searched wa s defined by the e xpected values in the Galactic plane 
given by the lCordes fc Laziol (2002) free electron density model of 
the Galaxy. Each data segment corresponds to one independent 
half of the original 35 minute integration. 



Parameter 


Value 


Integration time 


2x1050 s 


Sample interval 


1 ms 


Number of samples per data segment 


2 20 


Size of the data segments 


49 MB 


Acceleration search range 


±500 m s- 2 


Number of acceleration trials 


59 


DM search range 


0—1656 cm _3 pc 


Number of DM trials 


152 


Processing time (per data segment) 


~ 2.25 hr 


Total number of data segments 


2x34710 



iKramer et al.l 2006) and is a value that might be expected in 
compact NS-BH systems (see Section 6). Analyzing indepen- 
dent halves also provided better sensitivity to more compact 
binary systems (Pb > 3 hr) that would have been missed by 
acceleration searches of the full lengt h observation, b ecause 
of the effects of jerk (e.g. P b > 6 hr in lFaulkner et al.1 2004). 
The sampling interval of 1 ms was not optimal for millisec- 
ond pulsars, however it is expected that in most NS-BH 
systems the pulsar will be the second-born object and will 
not have undergone recycling (|Sipior et al.1 2004). 

Understanding the response of the data processing sys- 
tem to incorrect acceleration trials is important for validat- 
ing both genuinely accelerated and even solitary pulsar can- 
didates. The full effects of a constant line-of-sight acceler- 
ation on the integrated pulse S/N can be examined if we 
consider a simple Gaussian pulse convolved with the pulse 
broadening function given by Equation 2. In Figure [2] the 
bottom panel shows the signal advance or retard due to a 
number of constant accelerations within the 1050 s of the 
PMACCN search algorithm. The thin black lines represent 



1 




pulse broadening time (s) 

Figure 2. The bottom panel shows the signal advance or retard 
due to a number of constant accelerations. The thin black lines 
represent accelerations for which the data can be corrected us- 
ing the resampling method. The step size between acceleration 
trials is that given by Equation 3, for a sample interval of 1 ms 
and an integration time of 1050 s, as in the PMACCN search 
algorithm. In this case the simulated pulsar signal has an accel- 
eration of 78.3 m s -2 and falls in the middle of two acceleration 
trials (thick black line). The upper panel shows the effect of the 
pulse broadening on a simulated Gaussian pulse for a pulsar with 
a spin period of 50 ms and a pulse width of 10 per cent. The dot- 
ted line shows the original pulse profile, with no broadening; the 
thick black line shows the profile broadening caused by the accel- 
eration and the thin black line shows the "acceleration-corrected" 
profile for acceleration trials closest to the true acceleration. 



accelerations that can be compensated for using acceleration 
steps based on Equation 3. Here we plot the "worst case" 
scenario where a true signal from a pulsar, accelerated to 
78.3 m s -2 , shown by the thick black line, lies exactly in the 
middle of two acceleration trials. If resampling is performed 
with respect to the midpoint of the data then the maxi- 
mum pulse broadening is 4 ms. In the top panel we plot 
the effects of the pulse broadening on a simulated Gaussian 
pulse from a pulsar with a spin period of 50 ms and a pulse 
width of 10 per cent. The dotted line shows the original 
pulse profile with no acceleration effects. The thick black 
line shows the full effects of the pulse broadening caused by 
acceleration and the thin black line shows the 'acceleration- 
corrected' profile. The reduction in pulse amplitude, caused 
by broadening, depends on the intrinsic width of the pulse 
i.e. assuming an constant duty cycle, shorter period pulsars 
will experience increased broadening compared to longer pe- 
riod pulsars. In Figure [3] we show how the pulse S/N (re- 
lated to the effectiv e width and period of the pulse, see e.g. 
lLorimer fc Kramer! (2005) page 129) changes with incorrect 
acceleration trial for a number of pulse periods, ranging from 
two milliseconds to two seconds. In these examples the pulse 
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Figure 3. Relative pulse S/N variation with incorrect accelera- 
tion trial for a range of pulsar periods. In all cases the pulse duty 
cycle is five per cent, and the signal has a true acceleration of 
m s~ 2 . 
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Figure 4. Spectral S/N as a function of trial acceleration for four 
known relativistic binary pulsars using the PMACCN search al- 
gorithm. The horizontal dotted lines (at S/N = 7.4) show the 
minimum value of S/N considered not to be from the effects 
of statistical noise , based on the false-alarm probability (e.g. 
lLorimer fc Kramerl 2005). 



duty cycle is fixed at five per cent and the actual acceler- 
ation of the signal is m s~ 2 . Shorter period pulsars are 
smeared by incorrect acceleration trials more quickly than 
their longer period counterparts. As well as showing the im- 
portance of correcting for pulse broadening caused by con- 
stant accelerations, the curves displayed in Figure [3] provide 
a useful diagnostic for genuinely accelerated pulsar signals. 



3.1 Algorithm tests 

The PMACCN acceleration search algorithm has been 
tested on observations of a number of known relativistic bi- 
nary pulsars, and using d ata from s imulations of compact 
DNS and NS-BH systems (|Eatoughl 2009). Post facto tests 
of the algorithm, described in Section 16.31 have also been 
performed. The results from tests on real observations of rel- 
ativistic binary pulsars are presented in Figure[4] Each panel 
shows the spectral S/N, from the fluctuation spectrum, as 
a function of trial acceleration for 1050 s integrations with 
the PMPS observing system. PSR A in the Double Pulsar 
system is shown in panel (a) and is found at an accelera- 
tion of 122 m s~ 2 . Panel (b) shows PSR J1756-2251, de- 
tected at an acceleration of 52 m s~ 2 . In panel (c) the binary 
pulsar, B1913+16 is found with an acceleration of —26 m 
s~ 2 . Finally panel (d) shows the binary pulsar J1906+0746 
detected at —70 m s~ 2 . Standard non-accelerated pulsar 
searches would give spectral S/Ns roughly equal to the value 
in the 0ms -2 acceleration trial. In these tests the PMACCN 
acceleration search increases the spectral S/Ns by a factor 
of between at least two and five compared to a standard 
non-accelerated pulsar search. However, this improvement 
is dependent upon the orbital phase at which the pulsar was 
observed and is discussed in more detail in Section 6. 



4 THE PROCESSING 

It was estimated that processing the PMPS using the ac- 
celeration algorithm described in Section 3 would require an 
order of magnit ude more computational operations than the 
iFaulkner etail (2004) processing used. An opportunity arose 
to use the computing facilities of The University of Manch- 
ester, High Energy Physics (HEP) group. Their TIER2 com- 
puting facilitjjf] consists of 900 dual-core nodes that form 
part of a network of high-performance computers for the 
analysis of particle physics data from the Large Hadron Col- 
lider, known as GridPFfl GridPP was built by a collabora- 
tion of 19 UK universities, the Rutherford Appleton Lab- 
oratory and Conseil Europeen pour la Recherche Nucleaire 
(CERN). The UK's GridPP is connected to a much wider 
global collaboration of computing centres called the World- 
wide LHC Computing Griqj and in total provides the equiv- 
alent computing power of approximately 200, 000 CPUs. 
One of the primary aims of computational grids, such as 
GridPP, is to exploit and share the computational power 
available from disparate sets of computers or computer clus- 
ters that are connected by the Internet. Such 'grid' based 
distributed computing lends itself well to the highly-parallel 
nature of processing the data from pulsar surveys. 



4.1 Data processing overview 

Data processing was divided into two main stages. Firstly 
'local' processing, where the PMPS pointings were sepa- 
rated into 13 individual beams and each integration divided 
into two halves on the Jodrell Bank Centre for Astrophysics 



(JBCA) DCORE computer cluster. This was followed by 
'external' beam processing, where each half integration was 
acceleration-searched on the TIER2 facility. The flow of data 
and processing steps from acquisition at the telescope up 
to the inspection of pulsar candidates generated from the 
search is depicted in Figure [S] 

In summary, the PMPS data were stored in pointings 
of 13 independent beams on a RAID system at the JBCA. 
Each data file consists of 2 23 single-bit samples for each of 
the 96 3-MHz-wide frequency channels. The raw data were 
first separated into 13 individual beams using a dedicated 
FORTRAN program, SC_td. Each beam was then converted 
into multi-channel SIGPROC format, for use with the SIG- 
PROC pulsar processing suitjfl, and split into two 1050 s seg- 
ments. Each of the half-length integrations was stored on the 
JBCA RAID system, ready for processing on the TIER2 fa- 
cility. These initial processing tasks were performed on the 
DCORE computer cluster, to avoid excessive levels of data 
transfer to the TIER2 facility. 

The submission of processing jobs to the TIER2 facil- 
ity was controlled by the 'user interface' to GridPP: a linux 
desktop machine located at the JBCA. Processing jobs were 
submitted to the TIER2 'gatekeepers'; machines that main- 
tained a batch queue of jobs. Once a processing job started 
on a TIER2 worker node, a request was sent to a JAVA 
control program, running at the JBCA, that kept a queue 
of un-processed data files. When a data file became avail- 
able, it was copied to the TIER2 worker node. The data 
were then searched for periodic si gnals following th e stan- 
dard se^xcli_rjroredjrres_outlined in lManchester et al.l (2001) 
and lLorimer fc Kramerl (2005): Firstly, the 96 channels 
were dedispersed to four frequency sub-bands at each value 
of trial dispersion measure (DM). Just prior to dedisper- 
sion the harmful effects of radio frequency interference 
RFI) were red uced by application of the zero-DM filter 
Eatough et al.l 2009) . These sub-band data were stored for 
folding and optimisation of the candidates generated from 
the acceleration search. The subbands were then dedispersed 
(without re-application of the zero-DM filter) to a single 
channel timeseries at the same trial DM. Every four sequen- 
tial samples in this time series were then summed together, 
to give a 1 ms sample interval before performing the coherent 
acceleration search as outlined in Section 3. Features above 
a pre-defined threshold S/N (calculated from the expected 
false-alarm probability) in the resulting fluctuation spectra 
were stored in a candidate list along with the associated 
DM and acceleration. At this point 'reconstructed S/Ns' 
were also computed and storecQ. Harmonic summation of 
two, four, eight and 16 har monics was performed to provide 
sensitivity to narrow pulses jTaylor fc Hugueninl 1969). The 
threshold S/N was reduced in successive harmonic folds as 
the probability density function of spectral features evolves 
from a Rayleigh to a Gaussian distribution. The whole pro- 
cess was repeated for each trial value of DM and accelera- 
tion. 
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6 http:/ /sigproc. sourccforgc.net/ 

7 Reconstructed S/N is the S/N of a pulse profile gener- 
ated from the inverse Fourier transform of th e complex fre- 
quency components of a signal (e.g. Appendix C, iFaulknerl 2004, 
lLorimer fc KrameJ 2005). 
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Figure 5. Flow diagram summarizing the PMACCN data processing using the TIER2 facility. The survey pointings were stored on 
an archive RAID at the Jodrell Bank Centre for Astrophysics (JBCA). Each pointing was separated into 26 half integrations using the 
DCORE cluster and stored for processing. Jobs were submitted to the TIER2 facility via the 'user interface', a desktop machine at the 
JBCA. Once a processing job started on a worker node a request for data was made to a control program running on the JBCA RAID 
host computer. Available files were then copied and processed on the worker nodes. The candidates from the s earch were sent back to 
the DCORE cluster where they are viewed with either a graphical selection tool or an artificial neural network l|Eatough et alj 2010). 



After dispersion and acceleration searching was com- 
plete the candidate list was checked for the most significant 
signals (sorted by their reconstructed S /N) at each DM and 
acceleration. The DM and acceleration for which each sig- 
nal was strongest was noted and any harmonics of the signal 
were checked for and ignored if identified. Signals down to 
a threshold S/N ~ 9, from each of the 1050 s segments, 
were stored in a final candidate list. Candidates in this list 
were then 'optimised' in the time domain using the pul- 
SARHUNTER software packag^f). The appropriate sub-band 
data file was divided into 64 temporal sub-integrations after 
which searches over narrow ranges in period, DM and ac- 
celeration were performed to find the best integrated pulse 
S/N. The results of this search and the initial acceleration 
search were placed into 'candidate plots' for viewing. An ex- 
ample of one such candidate plot, covering the position of 
PSR J1141-6545, is visible in Figured Finally, all results 
were compressed with GZIP and returned to the JBCA for 
inspection. 



http: / / sourceforge.net/projccts/pulsarhuntcr/ 



4.2 Post-processing overview 

The total number of pulsar candidates produced by repro- 
cessing the PMPS with the PMACCN search algorithm 
was in the region of 16 million; a number far in ex cess of 
previous analyses. As first noted by iFaulkner et all (2004), 
such large numbers of pulsar candidates can, in practice, 
only be inspected with the aid of a graphical selection 
tool. For this purpose the graphical selection tool jreaper 
l|Keith et al.l 2009) was employed. Despite the advanced fil- 
tering techniques that can be performed with jreaper the 
number of candidates requiring visual inspection still re- 
mained large. As such, an automatic candidate inspection 
technique, utilising an artificial neural network (ANN), was 
developed to c omplement the manual searches performed 
with jreaper (|Eatough et al.l 2010). 

Promising pulsar candidates selected with either 
jreaper or by the ANN were placed into three categories for 
follow-up observation: class one, for candidates that had a 
large probability of being new pulsar discoveries, class three 
for candidates that had a smaller chance of being new pul- 
sars, and class two, for candidates that were in between the 
former categories. The majority of candidates viewed were 



File: PM0098_01761_21_OOl.phex.gz RA: 11:41 :05.8 Dec:-65:57:20.6 Gl:-64.15 Gb:-4.05 MJD:51465.85 
ObsFreq:1 51 6.5MHz Tobs:1 050.0s SourcelD:G4725480 Telescope:PARKES 
SpecSNR:62.4 ReconSNR:58.4 H-Fold: 1 6 
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Figure 6. Example candidate plot from the PMACCN algorithm generated for visual inspection of the quality of pulsar candidates and 
showing a detection the rclativistic binary pulsar, J1141— 6545. Starting clockwise from the bottom left the candidate plot displays: the 
integrated pulse profile, folded at the initial spectral discovery period and dispersion measure (DM); 64 temporal subintegrations of the 
observation, showing how the pulse varies with time; the spectral S/N as a function of a wide range of trial DMs (the DM-S/N curve); 
the spectral S/N as a function of trial acceleration (the acceleration-S/N curve: here the pulsar has been detected with an acceleration 
of 64.4 ms~ 2 ); the period-DM diagram, that shows how the S/N varies with small changes in the folding period and DM; the basic 
discovery parameters before and after time domain optimisation; stacked pulses across four frequency subbands, showing how the pulse 
varies with observing frequency, and finally the integrated pulse profile, folded at the optim um period and d ispersion measure. The dips, 
visible at the edges of the integrated pulse, are caused by application of the zero- DM filter ijEatough et al.l 2009). 



non-pulsar signals caused by either RFI or noise fluctua- 
tions. These signals were discounted. 



5 DISCOVERY AND TIMING OF 16 PULSARS 

Reprocessing of the PMPS with the PMACCN algorithm 
has yielded 16 new pulsar discoveries. Despite the improved 
acceleration search algorithm used in this analysis, no previ- 
ously unknown relativistic binary pulsars have been found. 
However, all known relativistic binary pulsars in the survey 
region have been re-detected with a higher significance than 
obtained in previous analyses. A number of promising binary 
candidates have been identified, but to date none have been 
re-detected despite follow-up observations (see Section 6.2). 

Not all of the pulsars found in this work have been de- 
tected here for the first time. During comparison of the can- 
didates generated from this search with those from the pre- 



vious [Rndkn^T^eFl^ (2004) analysis, it was noticed that a 
number of the previous class one candidates had never been 
re-observed. These oversights were undoubtedly due to the 
large number of genuine pulsar discoveries making tracking 
of the status of each candidate a difficult task. As part of this 
work, confirmation observations of these sources were per- 
formed after which they were added to timing programmes 
at either the Parkes or Jodrell Bank Observatories. All the 
pulsars have been timed for a period long enough to de- 
termine an accurate position and to derive the basic pulsar 
parameters. All timing measurements were performed with 
the PSRTIME softwar^. 

Table 2 presents the pulsar name, the J2000 right as- 
cension and declination from the timing solution, the Galac- 
tic coordinates, the S/N of the discovery observation af- 

9 http:/ /www .jb.man.ac.uk/~pulsar/observing/progs/psrtime. html 



ter the time-domain folding and optimisation, the candi- 
date selection tool by which the pulsar was discovered, 
the mean flux density at ~1.4 GHz from the survey de- 
tection, S1400, and the pulse width at 50 per cent of the 
peak of the mean pulse profile, W50. Uncertainties on the 
values presented, where relevant, are given in parentheses 
and refer to the last digit. Flux densities have been cor- 
rected for positional offsets in the survey detection. Sources 
marked with a * were not detected in this analysis but were 
found by reobservation of the unconfirm ed class one can- 
didates from the previous iFaulkner et~ai1 (2004) reprocess- 
ing. Sources marked with a * were detected both in this 
analysis and the previous analysis but only confirmed now. 
Sources with no markings have been detected here for the 
first time. In summary, five of the pulsars were not de- 
tected in this reprocessing, three were detected in both, 
and seven have been newly discovered. The discoveries of 
PSRs J1539-4828, J181 9-1717, and J1 835-0114 have al- 
ready been mentioned in lEatough et all (2009) and the dis- 
covery of PSR J1926+0737, using an ANN, is described in 
lEatough et~all (2010). 



Table 3 gives barycentric pulse periods, period deriva- 
tives, epoch of the period, the number of pulse time of arrival 
(TOA) measurements in the timing solution, the MJD range 
of the timing observations, the final rms residuals of the tim- 
ing model, and the dispersion measure. Table 4 gives the fol- 
lowing derived parameters of the 16 pulsars: characteristic 
age, t c , surface magnetic field, B 3 , rate of loss of rotational 
energy, E, the dispersion derived dis tance using the N E2001 
electron density distribution model (|Cordes fc Laziol 2002), 
the corresponding height above the Galactic plane, z, and fi- 
nally the inferred luminosity at 1400 MHz, Lmoo- High S/N 
mean pulse profiles of each of the 16 pulsars at 1374 MHz 
are presented in Figure [7] 



Excluding the pulsars that did not appear in this pro- 
cessing (i.e. sources marked with *), the remaining sample 
exhibits a wide range in flux densities. The discovery of PSR 
J1716— 4005, a source with a large flux density that was com- 
pletely undetected in previous analyses, has most likely been 
enabled by effective removal of RFI, by use of the zero-DM 
filter, and improved candidate selection tools. A number of 
the sources are below the limiting flux density of our search 
(5*1400 ~ 0.2 mjy for the centre of the central beam) which is 
estimated to be of order y2 times less sensitive than that in 
the original searches of the PMPS because of the use of half 
length, 1050 s integrations. The zero-DM filter may account 
for this as its use can improve the limiting survey sensitiv- 
ity by the effective removal of RFI (|Eatough et al.l 2009). In 
general the pulsars are consistent with the sample of pulsars 
found by the PMPS as a whole viz. low flux density but 
high luminosity because the pulsars are found at relatively 
large distances. PSR J1818— 1448 has a the largest period 
derivative of the pulsars discovered implying an age of just 
720 kyr and making it the youngest of the new sources. Un- 
fortunately, this age is too great to expect any supernova 
remnant association to be visible. 



5.1 Notable discoveries 

5.1.1 PSR J1835—0114, a binary millisecond pulsar 

Soon after the discovery of PSR J1835— 0114, it was noticed 
that the spin period was not constant but changing period- 
ically. As no sign of period variation was measured in the 
search output, i.e. over 1050 s integration, the behaviour 
was thought to be typical of that of a pulsar in a long or- 
bital period binary system. Continued timing observations 
at Jodrell Bank Observatory have shown that the pulsar is 
in a 6.7-day orbit, probably around a low mass companion. 
The measured binary parameters of this system are given 
in Table 5. The pulsar is in a circular orbit (e~lx 10 _ ), 
typical of pulsars in binary systems that have undergone a 
long period of accretion powered spin up (Phinney 1992). In 
non-relativistic systems estimation of the mass of the pul- 
sar and companion, m p and m c , can only be made through 
accurate measurements of the mass function given by 

r, • \ (m c sini) 3 4-7T 2 a 3 sin 3 i „„„„,,„„„„, 

f { m p ,m^) = ^-^L, = — — = O.OO2417O9Af , 

(4) 

where a is the semi-major axis of the orbit, i is the inclina- 
tion angle and the other terms have their usual meanings. 
In order to estimate the mass of the companion a typical 
neutron star mass of 1.35 Mq, from the known neutron 
star masses measured in relativistic binaries (Thorsett & 
Chakrabaty 1999), has been assumed. Random orbital axis 
orientation arguments allo w us to assume a med ian inclina- 
tion angle of 60°, see e.g. iLvne fc Smith! (|2006T ). page 127. 
With the observed value of the mass function, an iterative 
solution for m c in Equation [3] gives a most likely companion 
mass of 0.21 Mq; indicative of a white dwarf star. 

5.1.2 PSR J1808-1517, an intermittent pulsar 

PSR J1808— 1517 was a class one candidate that was de- 
tected in only one half of the original 35 minute observation. 
Analysis of the full integration, revealed the pulsed signal 
disappearing approximately half way through the observa- 
tion (see Figure [8] top panel). The candidate class, com- 
bined with observed parameters typical of normal pulsars 
(P ~ 544 ms, DM ~ 205 cm~ 3 pc, pulse duty cycle ~4 per 
cent) pointed to a clear pulsar detection. However, in ap- 
proximately seven hours of follow up observations performed 
on different dates at the Parkes observatory and regular ob- 
servations from March 13th to September 9th, 2008 at the 
Jodrell Bank observatory, no re-detection was made. Despite 
the expected difficulty in confirmation it was decided to con- 
tinue with observations at the Parkes observatory at approx- 
imately monthly intervals. In March 2010 a re-detection was 
made. To confirm the intermittent nature of this pulsar, re- 
observations were performed at the Effelsberg observatory. 
The bottom panel of Figure [5] shows temporal subintegra- 
tions over a 31 minute observation done at Effelsberg. Dur- 
ing this observation the pulsar exhibits nulls lasting the or- 
der of a few minutes, with a nulling fraction of ~ 50 per 
cent. Because of the earlier non-detections, we suggest the 
length of nulls vary and could be at least as long as a typi- 
cal observation time (> 30 minutes). Based on the Effelsberg 
observation the switching time scale between states is < 12 
s (< 23 pulse periods). 



Table 2. Positions, survey beam in which the pulsar was discovered, candidate selection tool used in the discovery, folded 
S/N, flux densities and pulse widths of the 16 pulsars discovered in this work. Positional uncertainties are derived from the 
one-sigma error of the timing model fit obtained using the PSRTIME software. See Section 5. for an explanation of the 
symbols. 
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Table 3. Timing parameters and dispersion measures of the 16 pulsars discovered in this work. Spin period uncertainties 
are derived from the one-sigma error of the timing model fit obtained using the PSRTIME software. See Section 5. for an 
explanation of the symbols. 
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This analysis of the PMPS, that has searched indepen- 
dent halves of the original 35 minute integrations, has aided 
in the detection of this pulsar. The addition of the second 
half of the original observation, when the pulsar was not visi- 
ble, reduces the integrated pulse S/N down to 8.5 and places 
the pulsar on the limit of previous detection thresholds. 



5.1.3 J1926+ 0737, a pulsar with an interpulse 

PSR J1926+0 737, that was found using an ANN 
l|Eatough et al.l 2010), exhibits an interpulse at a rotational 
phase ~ 180° from the main pulse. This characteristic im- 
plies that the angle between the magnetic axis and the 
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Figure 7. Mean pulse profiles over a phase of zero to one at 1.4 GHz for the 16 pulsars. The peak of each profile has been centered at 
a phase of 0.3. Each panel also displays the pulsar name, pulse period in seconds, and dispersion measure in cm - 3 pc. The effective time 
resolution of the profile, including the effects of interstellar dispersion, is indicated by the horizontal error bar below the text. 



rotational axis, a, is ~ 90° and that the pulsar must be 
an orthogonal rotator, where emission from both magnetic 
poles is visible. This pulsar adds to the population of 31 
normal pulsars ( with P > 20 ms) that exhibit interpulses 
(|Maciesiak et al.l 2011). The properties of pulsars with in- 
terpulses are of importance for understanding the evolution 
of pulsar beams, which has implications in the long term 
spin-down behaviour of pulsars and also in the number of 
pulsars likely to be det ected in future pulsar surveys (e.g. 
IWeltevrede fc Johnston! 2008, iMaciesiak et al.1 2011). 



6 DISCUSSION 



In this section we start with a brief comparison of our results 
with the available predicted number of relativistic binary 
pulsars in the PMPS (Section [67T}. We then outline some of 
the possible reasons for the non-detection of new relativistic 
binary systems (Section 16. 2[) . before describing further tests 
of the acceleration algorithm used here (Section l6.3|) . Lastly, 
Section [63] is a brief discussion of the promising pulsar can- 
didates detected with PMACCN, but which have never been 
confirmed in re-observations. 



Table 4. Deriv ed parameters of the 16 pulsars discovered in this work. Distances have been calculated based on the 
ICordes fc Laziol (2002) free electron density model of the Galaxy. See Section 5 for an explanation of the symbols. 
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Figure 8. Phase resolved pulse amplitude as a function of time 
for the original PMPS observation of PSR J1808-1517 (top panel) 
and a re-observation performed at the Effelsberg observatory 
(bottom panel). Subintcgrations arc plotted every ~ 30 s and 
~ 12 s for the PMPS and Effelsberg observations respectively. 



Table 5. Binary parameters of PSR J1835— 0114. Uncertainties 
on the values presented are given in parentheses and refer to the 
last digit. 



Parameter 



Value 



Orbital period (days) 6.6925427(4) 

Projected semimajor axis x = (a/c)sini (s) 4.654423(7) 
Orbital eccentricity e (xl0~ 5 ) 1.1(3) 
Longitude of periastron ui (°) -3(1) 
Epoch of periastron (MJD) 54540.7(3) 

Mass function (M ) 0.00241709(4) 



6.1 Comparison of results to binary population 
predictions 

To date, this acceleration search of the PMPS has not in- 
creased the number of relativistic binary systems discovered 
in the survey. Prior to the commencement of this work, and 
using a method b ased o n the observed sample of DNS sys- 
tems (|Kim et al.l 2003), iKalogera et al.l (2004) made a pre- 
diction that appro ximately anothe r 4 + 2 (error estimated 
from Figure 2. of IKalogera et all (2004)) compact binary 
pulsars (with properties similar to either PSRs B1913+16, 
B1534+12, and J0737-3039A/B) remained to be discovered 
in the PMPS data, and would be detecta ble with effective 
acceleration searches. The lFaulkner et al.l (2004) re-analysis 
of the PMPS proved that there were indeed binary sys- 
tems that had previously been missed because acceleration 
searches were not used, including one compact DNS system, 
PSR J1756-2251. PSR J1906+0746 was also detected with 
high significance in t his analysis bu t missed at the candi- 
date selection stage (|Lorimer et al.l 2006). Currently PSR 
J1906+0746 is believed to be a young pulsar with a heavy 



white dwarf companion, so does not belong to the category 
of systems described by Kalogera et al. (2004). The number 
of relativistic DN S systems discov ered in the PMPS is there- 
fore less than the iKalogera et aU (2004) prediction. If PSR 
J1906+0746 is in fact a DNS system, the current number 
of compact DNS in the PMPS would be two; still less than 
predicted. Whatever the nature of the companion, the youth 
of PSR J1906+0746 implies that the underlying popu lation 
of this type of system could be larger (|Lorimer et all 2006). 

Alternative Galactic population predictions, from 
binary synthesis methods, differ by orders of mag- 
nitude in the number of DNS and NS-BH systems 
produced because of large uncertainties in the model 
parameters used to describe key stages of binary evo- 
lution; such as the magnitude of supernova kicks and 
of com mon envelope evolution (se e e.g. 
1998, IPortegies Zwart fe Yungelsonl 1998, 
2001^ ISipior fc Sigurdssonl 2002, 
2002, 
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Belczvnski et all 2002, I Voss fc Taurisl 2002, 
Sipior et al.l 200 4. iLipunov et al.1 2005 . IPfahl et al.l 2005, 
Belczvnski et al.1 2010. lOslowski et al.1 2011). The popu- 



lation constraints provided by the binary pulsar search 
analysis performed here, and similar binary searches, can be 
used to improve and calibrate binary population synthesis 
codes. However, the result of a non-detection should be 
treated with caution due to a number of selection effects 
and stages at which pulsar candidates can easily be lost. 
Some of these points are discussed in the following section. 



6.2 Possible reasons for the non-detection of new 
relativistic binaries 

Although this search can help to constrain the numbers of 
new relativistic binary pulsar systems, our non-detection 
does not necessarily imply that there are no more of these 
systems to be found. A more likely scenario is that these 
pulsars have been missed because of any one of a number 
possible reasons, some of which are listed below. For exam- 
ple, such pulsar systems might be 

• below the limiting sensitivity of the PMPS. For this 
search of 1050 s segments, the best sensitivity of ~ 0.2 mjy 
is achieved in the centre of the central beam. 

• at an orbital phase where the acceleration was not con- 
stant (see Section |6. 3. 1[) . 

• outside the range of orbital accelerations searched (see 
Section EHO]) . 

• undetected because they are in a survey beam with ex- 
cessive RFI. Despite the use of the zero-DM filter, which 
reduces the amount of t he fluctuation s pectrum that is re- 
moved to virtually zero (|Eatough et al.1 2009), RFI can re- 
duce the sensitivity of t he survey. Unfortu nately it is difficult 
to quantify this effect ^Manchester et all 2001). 

• detected but not fo und during can didate selection (see 
e.g. iKeith et al.1 2009, lEatough et~ai1 2010). Pulsar sur- 
veys produce many millions of pulsar candidates which 
are typically inspected by eye and which can lead to 
promising candidates being overlooked. The missed dis- 
covery of PSR J1906+0746 in the PMPS illustrates this 
(|Lorimer et al.l 2006). Attempts in this work were made to 
automate this process using an ANN (see Section 4.2 and 
lEatough et~aH 2010). 



• already detected but cannot be confirmed with follow- 
up observations either because of intermittent emission (see 
e.g. Section I5.1.2|l or because the emission beam has pre- 
cessed out of the line-of-sight (see Section l6.4p . 

• undetected as they were not emitting during the sur- 
vey observation because they ar e a different class o f pulsar, 
e.g. Rotating Radio Transients dMcLaughlin et al.l 2006) or 
intermittent pulsars (|Kramer et al.l 2006) . 

• not in the region of sky covered by PMPS. See below 
for details of the new all-sky pulsar surveys. 

As already mentioned, some of these points are un-avoidable 
and others cannot be addressed until new more sensitive pul- 
sar surveys, such as th e all-sky Hig h Time Resolution Uni- 
verse (HTRU) surveys (|Keith et al.l 2010, Barr et al. in prep) 
are fully searched for binaries (|Nd 2011. Ng et al. in prep), 
or until surveys to be performed with the next generation 
of radio telesco pes, e.g. LOFAR, FAST, and the SKA (e.g. 
IStappers et al.l 2011, van Leeuwen & Stappers 2010, Smits 
et al., 2009a, b, Cordes et al., 2004) are completed. Some of 
the other points, relevant to this analysis, are now discussed 
in more detail. 



6.3 Additional acceleration search algorithm tests 

After th e pr evious re-proc essing of the PMPS 
(|Faulknerl 2004, iFaulkner et al.1 2004) it was understood 
that much of the sensitivity of the survey to relativistic 
binary pulsars was being lost because of two principal 
effects: firstly, the incoherent addition of spectra in the 
stack acceleration search and secondly, no compensation for 
the effects of jerk caused by large changes in orbital phase 
throughout the 35 minute observations. These problems 
have been partially addressed in this work by implementing 
a coherent acceleration search over a wide range of trial ac- 
celerations, and upon independent halves of the original 35 
minute integrations. To establish if any systematics in the 
PMACCN algorithm were responsible for the non-detection 
of any new relativistic binary pulsars we have performed 
a number of additional tests, which are described in this 
section. Because the primary goal of this work has been to 
find relativistic binary systems like the Double Pulsar (and 
those more extreme binary systems) we have chosen this 
system as our test case. 



6.3.1 Integration length analysed and orbital phase effects 

Since the acceleration search technique can only correct 
for the quadratic drift in pulse phase (and the correspond- 
ing linear drift in pulse frequency) caused by a constant 
acceleration, sensitivity would have been lost to any sys- 
tems that showed changes in acceleration during the PMPS 
observation. Without performing the higher order correc- 
tions necessary to remove the effects of jerk, the only 
method by which this can be mitigated in an accelera- 
tion search is to redu ce the integration length analysed 
l|johnston fc Kulkarnil 1991). In addition, the orbital phase 
at which the pulsar is observed will determine how well the 
constant acceleration assumption can model the data. To 
study both the effects of the integration length analysed (to 
find if the 1050 s integration time used in this work was op- 
timal for the detection of relativistic binary systems), and 
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Figure 9. The results from acceleration searches, using different integration lengths, across a 2.6 hour Parkes observation of the Double 
Pulsar system (panels (a), (b), and (c)) and the same data, but modified with a false sampling interval half that of the original to effect 
a 1.2 hour binary (panels (d), (e), and (f)). Each panel shows the spectral S/N (of PSR J0737— 3039A) as a function of trial acceleration 
for acceleration searches incremented in orbital phase by 100 seconds (panels (a), (b), and (c)) and 50 seconds (panels (d), (e), and (f)). 
The left-hand column of panels (panels (a) and (d)) show acceleration searches performed with an integration length of 671 seconds, the 
middle column (panels (b) and (e)) 1342 seconds, and the right-hand column (panels (c) and (f)) 2684 seconds. The black bar at the 
bottom of each panel indicates the length of the integration time used in each acceleration search. The maximum and average S/N peak 
are summarized in Table Al and are given here as follows. Panel (a): 62.6, 43.7, panel (b): 66.9, 48.3, panel (c): 62.8, 32.4, panel (d): 
66.6, 48.3, panel (e), 62.8, 32.4, panel (f): 36.4, 29.1. Contours, plotted on top of the colour scale, mark spectral S/Ns of 10, 30, and 60. 



to probe the effects of orbital phase at which the pulsar 
is observed, the following tests have been performed: Ac- 
celeration searches using different integration lengths, and 
separated by short intervals in orbital phase, have been per- 
formed on archival data, from the Parkes radio telescope, 
that covers one full orbit of the Double Pulsar system. Ad- 
ditionally, to create a more rigorous test, the same data have 
been modified with a false sampling interval half that of the 
original 80 /j,s. This has the effect of halving both the pulse 
and orbital period, and doubling the apparent line-of-sight 
acceleration to ~ 520 m s -2 . 

Three integration times, that are similar to small fac- 
tors of the original PMPS integration time, have been inves- 
tigated: "full-length" 2684 second segments, "half-length" 
1342 second segments, and "quarter-length" 671 second seg- 
mentj"!. Appropriate step sizes in acceleration trail for each 

10 The integration times investigated do not match exact factors 
of the original PMPS observation time because of the use of a 
different data sampling interval, and a technical requirement in 
the tests for a total number of samples equal to a power of two. 



integration length have been computed using Equation 3. To 
investigate the effects of orbital phase on the different length 
acceleration searches, the starting point of each search has 
been incremented by 100 seconds for the original data, and 
50 seconds for the modified data across the entire orbit. 



In Figure [9] the results of these tests are presented. 
Each panel shows the spectral S/N of PSR J0737-3039A 
as a function of trial acceleration for an acceleration search 
started at a time after MJD 52916.76 given in hours. Fig- 
ure [5] panels (a), (b) and (c) show the results from the orig- 
inal, un-modified data. Results from these tests are also 
summarized in Table 6 in addition to results from tests 
accounting for the effects of jerk (Section I6.3.3|) . Accelera- 
tion searches of the quarter and half-length segments (pan- 
els (a) and (b) respectively) are comparable in their sensi- 
tivity over all orbital phases although there is a small im- 
provement in sensitivity for the half-length segments, which 
show both a higher maximum and average peak spectral 



S/N (S/N pcak ) across the orbif^l. Importantly, the average 
S/N pcak in the half-length segments is larger only by a fac- 
tor of ~ 1.1 compared to the quarter-length segments; less 
than the expected increase of y/2 given by the radiometer 
equation. This can be explained by a reduction in sensitivity 
in the half-length segments because of the increased effects 
of jerk. This deleterious effect is confirmed by the fact that 
the maxima in S/N pcak occur at integrations that cover or- 
bital phases where the acceleration is closest to constant, 
viz. at the peaks and troughs. The corresponding minima 
in S/N peak occur at orbital phases where the acceleration is 
changing most rapidly. In the quarter-length analysis (panel 
(a)), where the acceleration can be considered constant at 
almost all orbital phases, the maxima in S/N peak occur more 
evenly across the orbit. 

As expected, the effectiveness of acceleration analyses 
using the full-length segments (panel (c)), is strongly de- 
pendent upon orbital phase. Although a maximum S/N peak , 
comparable to those in the shorter analyses, can be achieved, 
the average S/N peak is significantly reduced by a factor of 
~ 1.5 with respect to half-length analyses. 

Figure [9] panels (d), (e), and (f) show the results from 
acceleration searches of the data modified with a false sam- 
pling interval and using the same integration length seg- 
ments as previously. In this data, where the effective orbital 
period is now 1.2 hours, only acceleration searches that are 
performed on the short quarter- length segments (panel (d)) 
have sensitivity across all of the orbit. Both the values of 
the maximum and average S/N peak in this search (66.6 and 
48.3 respectively) and the occurrence of maxima in S/N peak 
at orbital phases where the acceleration is changing slowly 
(peaks and troughs) make this the direct analogue of the 
half-length analysis of the un-modified data (panel (b)), but 
with an acceleration range expanded by a factor of two. 
The same situation applies to panels (c) and (e). There- 
fore, an analysis based on 336 s integrations on the modified 
data should repeat the results of the 671 s analysis of the 
un-modified data (panel (a)). Tests using full-length inte- 
grations on the modified data produce the lowest values of 
maximum and average S/N peak from all the tests performed. 
Because of the relative brightness of PSR J0737— 3039A, de- 
tections can still be made at most orbital phases, however 
compared to the optimum acceleration search of the mod- 
ified data (671 s segments) the average S/N peak is reduced 
by 40 per cent. 

From these tests it is clear that the effectiveness of accel- 
eration searches for the detection of 'Double Pulsar like' sys- 
tems in long integration lengths (T b a > 2000 s) are highly 
dependent upon the orbital phase at which the pulsar is ob- 
served. With the integration length used in this work (1050 
s) the sensitivity to these systems lies somewhere in between 
the tests displayed in Figure [9] panels (a) and (b), suggest- 
ing this search should be sensitive to such systems at most 
orbital phases. For systems like the modified Double Pulsar 
system (Pb = 1.2 hrs, maximum acceleration ~ 520 m s~ 2 ) 
the possibility of detection in this work is strongly depen- 
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Figure 10. Maximum linc-of-sight orbital acceleration as a func- 
tion of orbital period for circular edge on binary systems. The 
neutron stars have a mass of 1.3 Mq and the two black hole com- 
panions have a mass of 8 Mq and 10 Mq. The solid horizontal 
line shows the magnitude of the maximum constant acceleration 
that this search was sensitive to of ±500 ms~ 2 . The diamond 
points show the maximum acceleration of PSR J0737— 3039A in 
the Double Pulsar system and the "modified" Double Pulsar sys- 
tem described in Section 6.1.1. 



dent on the orbital phase at which the pulsar was observed; 
the search is most effective when the pulsar is maximally 
accelerated. 



6.3.2 Acceleration search range 

Neglecting the effects of jerk, and only considering the mag- 
nitude of the orbital accelerations that might be present 
in even more compact DNS and NS-BH systems, it can be 
shown that the search performed in this work would not 
be capable of detecting these systems. Figure [101 shows the 
maximum line-of-sight acceleration as a function of orbital 
period, for circular edge on binary systems, as given by Ke- 
pler's laws. In all curves the neutron stars have a mass of 
1.3 Mq and the two black hole companions have a mass 
of 8 Mq and 10 Mq. The solid horizontal line shows the 
magnitude of the maximum constant acceleration that this 
search was sensitive to of ±500 ms~ 2 . The maximum ac- 
celerations present in the Double Pulsar system and in the 
'modified' Double Pulsar system, described in Section 6.1.1, 
are indicated by the diamond points. 

Using the PMACCN algorithm the most compact cir- 
cular NS-BH system that might be detectable would in be in 
an orbit of about three hours. To have an improved chance 
of finding the most extreme NS-BH systems (Pb > 2 hr) the 
acceleration search range should reach at least ±1000 m s~ 2 
(dotted horizontal line). Performing an acceleration search 
over this range, and with a computational cost no more than 
that used in the PMACCN algorithm, could be achieved by 
segmenting the PMPS data into four 525 s segments per 
beam. 



11 S/N peak refers to the best spectral S/N achieved in the acceler- 
ation search. The maximum S /N peak is the highest value found in 
the acceleration searches across the orbit. The average S/N pcak is 
an average of all acceleration searches performed across the orbit. 



6. 3. 3 Including jerk 

From an extension of the acceleration algorithm presented 
in Section 3 a step-size in jerk can be computed. The pulse 



Table 6. Summary table of results from acceleration and acceleration and jerk searches on the Double Pulsar system (J0737— 3039A) 
and the "modified" Double Pulsar system (J0737— 3039A*). S/N peak refers to the best spectral S/N achieved in the acceleration and 
acceleration and jerk search. The number of acceleration (or acceleration and jerk) trials at each DM searched, '"trials/DM, is quoted for 
a downsampled sampling interval used in each test of 1280 /is for J0737— 3039A and 640 /is for J0737— 3039A* . 
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broadening due to jerk, Tj e rk> can be written, 

7j„k(t) = Ji 3 /6c, (5) 

where j is the value of line-of-sight jerk, typically given in 
units of cms -3 , t is time and c is the speed o f light (see 
for example Equation 4 in I Johnston fc Kulkarnil 1991). Fol- 
lowing the principles outlined in Section 3 the maximum 
acceptable value of pulse broadening due to jerk at either 
end of the integration should be 8r samp for at least 50 per 
cent of the pulses to be smeared by less than one time sam- 
ple. A maximum pulse broadening time of 8r B am P allows a 
step of 16r S am P between jerk trials. Setting Equation [5] equal 
to this value and once again letting t = T b s /2, the step-size 
in jerk, 5j can be written as, 

Sj = 768cT sa mp/Tobs- (6) 

To find the sensitivity to be gained by correcting for jerk 
we have re-performed the tests described in Section [6.3.1l on 
the un-modified Double Pulsar data (Figure[9]panels (a),(b), 
and (c)). Step sizes have been computed with Equations 3 
and [6] and a range in jerk of ±20 cms -3 (c.f. jerk in Double 
Pulsar system ~ 18 cms -3 ) has been searched. Again, the 
maximum and average values of S/N peak across the orbit for 
the different integration lengths are summarized in Table 6. 
Tests on the 'modified' Double Pulsar system have not been 
performed due to limitations in available computing power. 
For the the shortest quarter-length analysis (671 s) the re- 
sults are the same as the acceleration tests as no jerk trials 
are required within the range we have chosen. For the half- 
length segments there is no change in the maximum S/N peak 
achievable but the average S/N peak has increased by seven 
per cent. The biggest improvement from the application of 
jerk corrections is in the full-length (2684 s) segments where 
both the best maximum and average S/N poak from all of 
the searches are achieved. These results can also be seen in 
Figure [TT1 where the cumulative value of S/N pcak of PSR 
J0737— 3039A from the acceleration and acceleration and 
jerk searches across the Double Pulsar orbit are plotted. 

From Table 6 and Figure [11] it can also be seen that the 
improvement in average S/N pcak across the orbit between 
acceleration and jerk searches of the half-length and full- 
length integrations is small (factor of 1.04). Because there 
are a factor of ~ 21 times more trials to perform the accel- 
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Figure 11. The cumulative S/N pcak for detections of PSR 
J0737— 3039A from the acceleration and acceleration and jerk 
searches across one orbit of the Double Pulsar system. 



eration and jerk analyses on the full-length segments com- 
pared to the half-length segments, there appears to be large 
computational expense with little gain in sensitivity. These 
lower than expected gains in sensitivity by including jerk 
trials are most likely due to poor approximation of the line- 
of-sight orbital motion over the length of the integration. 
To combat these effects a third derivative of velocity might 
need to be included. 

6.4 Pulsar candidate non-confirmations 

Although no relativistic binary systems have been discov- 
ered the acceleration search has produced a number of 
high-quality candidates that have not been confirmed, de- 
spite multiple re-observations at either the Parkes or Jodrell 
Bank observatories. Figure [121 shows an example of one such 
source, J1431— 6042, a highly accelerated class three can- 
didate. This can didate has be en followed up with a 'grid' 
observation (e.g. iMorris et al.l 2002) and around four hours 
of re-observations at Parkes. Because the orbital phase and 
hence the orbital acceleration are likely to be different upon 
re-observation all the data have been searched in accelera- 
tion on the TIER2 facility with the PMACCN algorithm, 
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Figure 12. The class three relativistic binary pulsar candidate, J1431— 6042. Starting clockwise from the bottom left the candidate 
plot displays: the integrated pulse profile, folded at the initial spectral discovery period and DM; 64 temporal subintegrations of the 
observation, showing how the pulse varies with time; the spectral S/N as a function of a wide range of trail DMs; the spectral S/N as a 
function of trial acceleration (here the pulsar candidate has been detected with an acceleration of —255 ms -2 ); the period-DM diagram, 
that shows how the S/N varies with small changes in the folding period and DM; the basic discovery parameters before and after time 
domain optimisation; stacked pulses across 32 frequency subbands, showing how the pulse varies with observing frequency, and finally 
the integrated pulse profile, folded at the optimum period and DM. 



and independently with the presto software suit43 over 
a wider range of accelerations and including corrections for 
jerk. To date no conclusive re-detection of this candidate has 
been made. 

In total there are 12 other candidates for which a similar 
situation applies. The majority come from our class two and 
class three candidate lists. Typically these candidates appear 
accelerated, have a S/N near our detection threshold (S/N ~ 
9), and are often only found in one of the half integrations. 
Many of these detections might be attributed to strong 'de- 
accelerated' RFI signals. Strong periodic RFI can appear 
in the higher DM trials and sometimes even swept in radio 
frequency mimicking a dispersed signal. In addition, RFI 
signals are typically unstable causing them to drift through 
a number of spectral bins in the fluctuation spectrum. 

If any of these candidates are in fact true celestial 
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sources the systems may have undergone geodetic pre- 
cession, where the spin axis of the pulsar precesses with 
respect to the total angular momentum of the system, 
which is dominated by the orbital angular momentum 
jDamour fc Ruffinil 19741). In a manner similar to that shown 
for PSR B1913+16 bv lKramerl (jl99&t ). the effect may have 
shifted the orientation of emission beam away from the 
Earth over the ~ 10 years since the original PMPS obser- 
vation. This is the current expla nation for why the relativis- 
tic binary pulsar, J1141— 6545 (jKaspi et al ] 120001 ) was not 
detec ted in an earlier 7 cm survey that covered the posi- 
tion ([Manchester et all 2010). We can estimate the magni- 
tude of geodetic precession for candidates that are close to 
the orbital period limit of our search. Assuming a minimum 
detectable orbital period of ~ 3 hrs (using the PMACCN 
search algorithm) gives a precessional rate of 3.5° yr _1 for 
circular DNS systems (neutron star mass of 1.3 Mq) and 
14.7° yr _1 for circular NS-BH systems (black hole mass of 
10.0 Mq). Although the alignment between the spin axis 



of the pulsar and orbital angular momentum axis is impor- 
tant in determining by how much the emission beam will 
move, it is possible that relativistic systems detected in the 
PMPS survey observations may have precessed out of the 
line-of-sight in the intervening years. 

At the time of writing, observations to detect un- 
confirmed candidates from this search are still performed 
when telescope time becomes available. 



7 SUMMARY & FUTURE WORK 

A full re-analysis of the PMPS has been performed using 
coherent acceleration searches on independent halves of the 
original 35 minute integrations. 16 pulsars have been discov- 
ered, including a binary millisecond pulsar, an intermittent 
pulsar, and a pulsar with an interpulse. No new relativistic 
binary pulsars have been discovered. 

Searches for relativistic b inar ies are already underway 
in the HTRU pulsar surveys (|Ng| 2011, Ng et al. in prep). 
The low-latitude portion of the survey will cover a thin strip 
(|6| < 3.5°) along the Galactic plane with integration times 
of over 70 minutes. For these long integration times acceler- 
ation searches will not be sensitive to the most relativistic 
binary pulsars; higher order effects will need to be taken 
into account (see Section 16 . 3 . 3 p . Unfortunately, it can be 
shown that searches in these higher order effects are still 
computationally prohibitive for long integrations. Assuming 
a sampling interval of 128 /xs and using step-sizes in accel- 
eration and jerk given by Equations [3] and [6] respectively, 
corrections for the effects of both acceleration and jerk in 
a 'Double Pulsar like' system would require of the order of 
~ 10 7 trials per DM trial for a 70 minute observation (c.f. 59 
acceleration trials per DM as in this work). For the analysis 
of the low-latitude HTRU data, either incoherent methods 
such as stack searches (Eatough et al. in prep), or coherent 
acceleration searches of different length but shorter segments 
of the original 70 minute observations should be performed 
(Ng et al. in prep). 

Because of the limitations of the acceleration search 
performed in this work, a full re-processing of the PMPS 
has been conducted utilizing the Einstein@Home net- 
worlT^I. This search applies matched filters to coherently 
de-modulate the data based on a number of Keplerian or- 
bital templates, which allows the full integration length to 
be searched (Knispel et al. in prep). 

The discovery and timing of highly relativistic binary 
pulsars is on e of the key science projects of the SKA 
(|Cordes et al.l 2004). As well as offering highly reward- 
ing discoveries, current searches for such systems, both in 
archival data and in data from the latest pulsar surveys, 
will provide considerable information on the binary param- 
eter space still to be searched in, and on the strategy of 
future pulsar surveys. 
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